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There is a single negative mode in the spectrum of small perturbations about the tunneling solutions describing 
a metastable vacuum decay in flat spacetime. This mode is needed for consistent description of decay processes. 
When gravity is included the situation is more complicated. There are different answers in different reduction 
schemes. An approach based on elimination of scalar field perturbations shows no negative mode, whereas the 
recent approach based on elimination of gravitational perturbations indicates presence of a negative mode. In 
this contribution we analyse and compare the present approaches to the negative mode problem in false vacuum 
decay with gravity. 
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1. INTRODUCTION 

The first order phase transitions might play an 
important role in the Early Universe They 
proceed via nucleation of the bubbles of true vac- 
uum in metastable (false) one and subsequent 
growth of the bubbles. The false vacuum decay 
is usually discussed in frame of self-interacting 
scalar field theory. This process is described by 
the O (4)— symmetric bounce solution ||,^ of the 
Euclidean equations of motion. Value of the Eu- 
clidean action at the bounce gives leading ex- 
ponential factor in decay rate. The perturba- 
tions about the bounce solution define the one 
loop corrections to the bubble nucleation rate and 
determine quantum state of materialized bub- 
ble It is remarkable that in the spectrum 
of small perturbations about the bounce in flat 
spacetime there is exactly one negative mode 
This mode is responsible for making correction to 
the ground-state energy imaginary, thus justify- 
ing decay interpretation. 

When gravity is included the model contains 
gauge degrees of freedom and to extract informa- 
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tion about the spectrum one should find proper 
reduction to physical variables. In the Euclidean 
gravity the presence of the gauge degrees of free- 
dom manifests itself in the well known confor- 
mal factor problem. In pure gravity it is pos- 
sible to cure this problem either at the level of 
prescription Q, or via careful gauge fixing 
For scalar perturbations in the theory of a scalar 
field coupled to gravity it is the only spatially 
homogeneous modes, which suffer from the con- 
formal factor problem (see e.g. [Q). On the 
other hand the negative mode is expected exactly 
in this problematic sector of scalar homogeneous 
perturbations, which consists of coupled scalar 
field and metric perturbations. The task is to 
reduce this coupled system to a system with a sin- 
gle physical variable. Note that there is no prob- 
lem in elimination of unphysical variables in per- 
turbations about the flat Friedmann-Robertson- 
Walker (FRW) type model with scalar field [||. 
Problem arises for perturbations in the theory of 
scalar field in closed FRW type Universe. 

There are two possible ways of reduction of cou- 
pled system of perturbations: either to eliminate 
metric perturbations or perturbations of scalar 
field. The approach based on elimination of scalar 
field perturbations shows no negative mode in 
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the spectrum of remining variable ||^- [l^,|7 |. On 
the other hand the recent investigation |13| indi- 
cates presence of a single negative mode in the 
approach based on elimination of gravitational 
perturbations. This indicates on lack of full un- 
derstanding of the role of gauge-fixing procedure 
in description of false vacuum decay with gravity 
and need for further investigation of this subject. 

The aim of present contribution is to discuss 
and compare existing approaches to the negative 
mode problem. The rest of the article is orga- 
nized as follows. In the next section we recall 
main formulas relevant to false vacuum decay in 
flat spacetime. In Sec. 3 we discuss inclusion 
of gravity, remind some important classical solu- 
tions of Euclidean equations of motion and derive 
quadratic action for scalar perturbations about 
0(4)— symmetric background solutions. We will 
restrict ourselves to the most problematic sector 
of 0(4)— symmetric perturbations. In Sec. 4 we 
analyse three different reduction schemes. Sec. 5 
contains the concluding remarks. 




Figure 1. Schematic view of a scalar field poten- 
tial V{(l)) with a metastable vacuum (jj = (jj^, true 



vacuum 



= f/)_ and local maximum <h — 



2. FALSE VACUUM DECAY IN FLAT 
SPACETIME 

Let's consider self-interacting scalar field (f) in 
flat spacetime and assume that the scalar field 
potential T^((/>) has shape shown in Fig.l, with 
local minimum (false vacuum) at = and 
absolute minimum (true vacuum) at = 0_ . The 
false vacuum decay is described by the Euclidean 
path integral 
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-HT 



\<f>+ >-- 
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where T is a large positive number and H is the 
scalar field Hamiltonian 



H{tt,, 



^ + i9„0(x)a"</.(x) 



+V{(j))](fx, 1,2,3. 
Integrating Eq.(P over tt^ one gets 

-Se jj. 
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where Se is the usual Euclidean action 

Se^ f [ld^<l>{x)d'^^{x) + V{^)]d^xdT , (4) 



with n ~ 0,1, 2, 3. In the limit T oo the l.h.s. 
of Eq.(||) contains information about the lowest 
energy eigenvalue and wave function. The r.h.s. 
can be calculated in the WKB approximation. 
The energy Eq of the lowest state gets correction 
due to tunneling phenomenon 



E ^ Eo- Ke^^ 
with 
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where (pbounce is the 0(4)— symmetric euclidean 
bounce solution with the lowest action, which in- 
terpolates between true and false vacua. The fac- 
tor K is obtained by the Gaussian integration of 
exponent of quadratic action of small perturba- 
tions about the bounce 



Hop- 
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4^ \det[-d'^ + V"{(t)+)] 



-1/2 



(7) 
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The operator -9^ + V"{(j)) has zero Diodes and 
det' means that the determinant computed with 
the zero eigenvalues omitted. It turns out that 
there is exactly one mode with the negative eigen- 
value in the spectrum of small perturbations 
about the bounce solution. This makes K and 
hence the energy shift in the Eq.(||) purely imag- 
inary and supports the false vacuum decay in- 
terpretation. Decay rate per unit volume V, per 
unit time is given by 



r 

V 



\K\e- 



(8) 



Note that the functional integral Eq. is written 
for unconstrained (physical) degrees of freedom. 
If one has m first class constraints 



Ca{(t>,TT) W , a = 1, ...,TO , 



(9) 



then according to general recipy [|4j one should 
choose some gauge fixing conditions 



X;3((/>,7r) = , ^ = 1, ...,m , 



(10) 



and modify the functional integral, Eq.(|l|), as fol- 
lows 

J^'^^^-^^^'^A 5{C^) 5{x^)dc^diT^ , (11) 

with the Faddeev-Popov determinant 

A = det{C„,x''}^0. (12) 

3. INCLUSION OF GRAVITY 

The Euclidean action of system composed of 
scalar field minimally coupled to gravity is 



S 



-^i? + iv^0v^0 + y(0) 



(13) 



where n = SttG is the reduced Newton's constant. 

Since we are interested in solutions having 
0(4)— symmetry, the metric is parametrised as 
follows 



da^ + a^{a)^ijdx^dx\ cf) 



(14) 



where 7^ is the three-dimensional metric on the 
constant curvature space sections. The field equa- 
tions are 



S+y(r + ^('^))=0, 



(15) 
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where a dot denotes a derivative with respect to 
proper time a and /C is the curvature parame- 
ter, which has the values 1, 0, -1 for closed, flat 
and open universes respectively. In what follows 
we are interested in closed case, but we do not 
specify the value of K. now in order to see how it 
enters in the final expressions. 

3.1. Euclidean solutions 

There are few celebrated solutions of these 
equations, which play an important role in Eu- 
clidean quantum gravity. 

1. The Hawking-Moss solution is a 4-sphere 
jl^ ] corresponding to scalar field sitting on the 
top of the potential barrier 



(18) 



with Htop — \J KV^('/'top)/3. Note that the Eu- 
clidean time a varies in finite interval a = (0, cr/). 

2. The Coleman-De Luccia bounce is a de- 
formed 4-sphere ||^. It starts with some (j) = (po 
at (T = close to 0_, stops at a = erf close to 
and obeys the regularity conditions 



a(0) 0(0) = 0, a(cr/) = 0(cr/) = 



(19) 



3. Relaxing the regularity conditions at the 
second zero of a one gets the Hawking- Turok sin- 
gular instanton [ p6[ . It was introduced to describe 
creation of an open inflationary Universe. 

Our main goal is to investigate spectrum of 
small perturbations about the Euclidean solu- 
tions and, in particular, to determine presence 
of a negative mode. 

3.2. Quadratic action of scalar 0(4)— sym- 
metric perturbations 

Investigation of perturbations would be con- 
venient to perform in conformal frame. We 
expand the metric and the scalar field over a 
0(4)— symmetric background 

ds^ = a(T)2 [(1 + 2A{T))dT^ 
+ -/^j{l-2^{T))dx'dx^] , 
= ^(T) + $(r), (20) 
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where r is the conformal time, a and <p are the 
background field values and A, ^ and $ are small 
perturbations. 

The background field equations in the confor- 
mal time are 

-H' -IC = ^^'\ (21) 
2H' + H^-IC = -'^(>fi'^ + 2a^V{(p)), (22) 

if" + 2nip' -a'^^ = 0, (23) 

OLp 

where a prime denotes a derivative with respect 
to T and H := a' /a. 

Expanding the total action, keeping terms of 
the second order in perturbations and using the 
background equations, we find 



(24) 



where S^^^ is the action for the background solu- 
tion and S^"^^ is quadratic in perturbations with 
the Lagrangian for scalar O (4)— symmetric per- 
turbations 

2k '■ 

+ K($'2+a2i!^$2^g 

dipdip 

At/ 

-(2k^'$' - 2Ka2 — $ -1- 127^*' + 12/C*)A 



6ip 



-2{H' + 2H^ +IC)A^] . (25) 



Note that the variation with respect to A gives 
the first order (constraint) equation 



AT/ 

2Kifi'^' - 2/40^ — $ + 12H^' + 12/C^' 

dip 

+4:{n' + 2n^ + jc)A^o . 



(26) 



4. THREE APPROACHES TO THE RE- 
DUCTION 

To obtain the unconstrained system corre- 
sponding to the degenerate Lagrangian (|2^) we 
will follow the conventional Dirac formulation of 
generalized Hamiltonian dynamics Cal- 
culating the canonical momenta 

(5(")/: /7 / , K , 



^/7(cI>'-^'A), 



0, 



(27) 



we find the primary constraint Ci := Ha = 0. 

Thus the evolution is governed by the total 
Hamiltonian 



Ht^Hc + ui{t)Ci 



(28) 



with arbitrary function ui(t) and canonical 
Hamiltonian 



He 



1 



12a2/7 



2a2^ 
2 \ Sip5ip 



imp 



— Kip ■ 
2 ^ 



dip K 



(29) 



Conservation of primary constraint gives the sec- 
ondary constraint 



.dV 

Sip 



(30) 



The primary and secondary constraints are first 
class one and there are no ternary constraints. 

Under the infinitesimal shift r ^ r + A these 
constraints generate the gauge transformations 



6^ = -n\ , 

(5$ = ip'X , (5H$ 
6 A = X' +HX. 



6U, = Z!^!^x , 



ia^^iip" -ip'n)X 



(31) 



The existence of constraints in the system as 
usually means the presence of unphysical degrees 
of freedom and that the whole system has no 
unique dynamics. To find the physical variables 
with unique evolution one should fix the gauge 
and solve the constraints. There are two pos- 
sible general strategies: either to eliminate the 
gravitational degrees of freedom (H$ and ^E") or 
perturbations of scalar field (H^ and <&). 
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4.1. The I approach 

In the first investigation on this subject the La- 
grangian approach was used fl^ . The gauge was 
fixed by the condition 



* = 



(32) 



Using this gauge condition and ehminating A 
with the help of the constraint equation Eq.(26) 
we obtain the unconstrained quadratic action in 
the form 



S 



(2) 
LRT 



'^Qlrt 

— — + Q 

6 Sifi 



a2 ^ 



3 



Sip 



LRT 



fv_\ 

S(pS(p J 



with 



hRT in' + 2n^ + /c)/3 = 



drd^a 



Kip 
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(33) 



(34) 



We see that the action Eq. (^3|) has correct over- 
all sign if Qlrt > 0. As it was shown in 
under certain choice of parameters the bounce 
can develop a region with negative Qlrt- The 
main conclusion in JTot was that the region with 
Q LRT < leads to catastrophic particles creation 
and is pathological. 

However the gauge fixing Eq. (^) has problems 
where a' = 0. This can be seen in the Hamilto- 
nian formalism, while the Faddeev-Popov deter- 
minant is proportional to a' . 

4.2. The II approach 

Eliminating 11$ using the constraint Eq.(30) 
and working in the gauge invariant variables 

* = 'I'H $, 



n 



+ ^$ 

Kip' 



(35) 



one obtains the reduced hamiltonian 

-.2 



H* = -- 



K 
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K 



K 



in 



ip'- \ K 2a^^ 
Performing the canonical transformation 



(36) 



IKLO 



n 



3ICa^^/jip' 



iJCa^^yp'^ 2i 
= q - 

2n Kip 



(37) 



and solving for the momenta p one obtains [|10| 
the quadratic part of the Euclidean action 



^(2) ^ (1 - 4/C) 



Uq^ 



y^d^xdr , (38) 



with a potential U depending on the background 
fields 



K ,2 



4/C 



(39) 



and q — aq. 

We see that quadratic action for the homo- 
geneous harmonic has "wrong" overall sign. To 
overcome this problem it was suggested |9| that 
analytic continuation 



'iq 



(40) 



is performed while integrating over this mode 
(compare flpot). 

The equation for the mode functions, which 
diagonalize the action (^8|), has form of the 
Schrodinger equation 



dT^ 



q + Uq^ Eq 



(41) 
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Note that the singularities — , which appear 

ip' 

in this approach do not allow one to investi- 
gate straightforwardly the perturbations about 
the Hawking-Moss solution (with ip = Ptop = 
const). 

The spectrum about the Coleman-De Luccia 
bounce was investigated and no negative mode 
theorem was proven [ pT| , p^ . 

The Hawking- Turok instanton was investigated 
and it was found that for monotonous scalar 
field potentials there is no negative mode in the 
spectrum. On the other hand it was found that 
the negative mode is present in "exotic" cases 
when the Hawking- Turok instanton "overshoots" 
the Coleman-De Luccia bounce. This means that 
in this case the Hawking- Turok instanton is un- 
stable and has larger action then the correspond- 
ing Coleman-De Luccia bounce. 

Note that the formulation in terms of the gauge 
invariant variables Eq.(pq) corresponds to the 



6 



gauge choice xgmst := $ = (compare [[lO[). 
Hence the Faddeev-Popov determinant is ^p' and 
this approach is singular for certain configura- 
tions. 

4.3. The III approach 

Having in mind that i) spherically symmetric 
gravity has no propagating degrees of freedom 
and ii) in the limit k — *■ we should recover scalar 
field theory in flat space, we prefer an approach 
based on elimination of gravitational perturba- 
tions At the same time from the Eq.(|30|) 
we see that solving constraint we divide either to 
If' or to Ti or to a^. Since ip' and a' might have 
zero(s) by solving constraints we might introduce 
extra singularities in quadratic action. For closed 
universe the scale factor a is vanishing for some 
r = T±. Since we are interested in perturbations 
which vanish for t — > rj-, the coordinate singu- 
larity a(r±) — looks most harmless. 

We use the following gauge fixing conditions 



Xi 



0, 



X2 



:H, 



, (42) 



which obey the following canonical algebra: 

{C„Q} = {x.,x,} = 0, = (43) 

As a next step we consider all constraints in a 
strong sense and define the physical Hamiltonian 
as 



H* := Hr 



(44) 



For the physical Hamiltonian wc find 



H* = 



+ 



1 / ^^ i2\ 

-f (a^ 3^'7^)H*$ 
o/C Of) 



2 ^'Sifi'nf 



,fv_ 



> 12 



(45) 



which corresponds to the following unconstrained 
quadratic action for one physical dynamical de- 
gree of freedom 



2Q 



with 
Q:= 



1 



6/C 



(46) 



(47) 



In what follows we will consider the background 
configurations for which the factor Q is positive 
definite. 

hitroducing a new variable q = aj^fQ $ and 
integrating by parts we obtain the unconstrained 
quadratic action in the form 



2 



(<Z'2 + W{a{r), p{T)\q^) dr^d^x ,(48) 



with frequency W whose time dependence is de- 
termined by the background solutions 



W\a{T),p{T)\ = 



Q 6ip5ip 2)CQ'^ Sip 



/CQ2 Sp 



Q 



/C - 6/C - 3/cg 



(49) 



Note that the obtained action allows us to con- 
sider the slowly varying scalar field, p' 0, and 
vanishing gravity, k — > 0, limits. 

The equation for the mode functions, which 
diagonalize the action Eq.(|4^), has form of the 
Schrodinger equation 



-j^q + W[a{T),p{T)]q 



Eq 



(50) 



with the potential Eq.(|49|). Having this equation, 
one can investigate the negative mode problem 
for concrete cases numerically. We investigated 
the case with the scalar field potential 



V{p>)^—{p>\p>-vf+Bp^) 



(51) 



This potential for = 2, B = 0.12 and v = 0.5 
has local maximum at ptop = 0.31250, and local 
minimum (false vacuum) at (p false = 0.3571428. 



1 units are used. 



7 



For this potential there exists the Coleman-De 
Luccia bounce solution with ipo = 0.1123579. 
The quantum mechanical Schrodinger equation 
Eq.(50) was solved numerically for this potential 
and exactly one bound state was found. Note that 
the factor Q was positive definite for this case. 

Within this approach it is problematic to inves- 
tigate the Hawking- Turok instanton, while scalar 
field (fi runs away and eventually the factor Q be- 
comes negative. 

We found that the Schrodinger equation 



Eq.(50) for the Hawking-Moss background solu- 
tion has six boundstates. This number is in per- 
fect agreement with the analytic formula for the 
eigenvalues 

A„ = n{n + 3)nlp + V"{iptop), n = 0, 1, 2... (52) 



while for our example 



top) 



''-top 



40.96 . 



(53) 



Presence of many negative modes about the 
Hawking-Moss solution means that the corre- 
sponding Coleman-De Luccia bounce for the 
given potential and set of parameters gives dom- 
inant contribution to tunneling ||2^,^ . 

Note that this approach based on gauge fix- 
ing Eq. ( p2[ ) can be equally well formulated in the 
gauge invariant variables 



* = $ 



:n>. 



n* = n$ + — ((^"-(^'7i)n>p . (54) 

6/C 



indicator for understanding to which type each 
solution belongs: instanton (mixing) or bounce 
(decay). 

The approach based on elimination of gravita- 
tional degrees of freedom shows no negative mode 
in the spectrum. It was suggested j9| that the i 
which comes from the conformal rotation Eq. ( |4C| ) 
plays the same role as a negative mode. This ar- 
gument has a weak point. The conformal rotation 
Eq.(|40|) is the same for all euclidean background 
solutions, while the imaginary shift in energy is 
peculiarity of an unstable state. Moreover, for the 
Giddings-Strominger wormhole a negative mode 
was found on top of conformal rotation [ pO{ . 

The approach based on elimination of gravita- 
tional degrees of freedom indicates the presence of 
a negative mode. It has smooth k — > limit and 
can be straightforwardly used for the Hawking- 
Moss solution (with (p = 0). At the same time 
there are still many open questions. A detailed 
numerical study of different cases is needed. An- 
other issue which needs further investigation is an 
understanding of the role of configurations with 
the negative factor Q: is it real physical instabil- 
ity or it is the breakdown of the reduction scheme. 

I am deeply thankful to the organizers of QG99 
for splendid hospitality and for creating charming 
atmosphere for discussions during the meeting. 

Note added: After this contribution was essen- 
tially completed, further progress in investigation 
of negative mode problem was made. The results 
are summarized in the revised version of |l3|. 



5. CONCLUDING REMARKS 
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